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Abstract. We present measurements of CP violation parameter <fii//3 in B° decays that are domi- 
nated by b — > ccs, b — ► sqq and b — > dqq transitions. The results are based on a large sample of BB 
pairs recorded at the T(4S) resonance with the Belle detector at the KEKB energy-asymmetric 
e + e~ collider. CP violation parameters for each decay mode are obtained from the asymmetries in 
the distributions of the proper-time intervals between the reconstructed B and the accompanying 
B meson. 
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1 Introduction 

The Standard Model (SM) describes CP violation in 
B° meson decays using the complex phase of the 3x3 
Cabbibo-Kobayashi-Masukawa (CKM) mixing matrix 
H]. In the decay chain T(4S) -> B°B° - fcpftag, 
where one of the B mesons decays at time tcp to a final 
state fcp and the other decays at time t tag to_a final 
state ftag that distinguishes between B° and B°, the 
decay rate has time dependence [2! given by V(At) = 

e -\M/r B o 

— [l + q[S f sm(Am d At) + A f cos(Am d At)}] . 

Here Sf and Af are CP-violation parameters, t b <> is 
the B° lifetime, Am<i is the mass difference between 
the two B° mass eigenstates, At = tcp ~ itag> and 
the 6-flavor charge is q = +1(— 1) when the tagging B 
meson is a B°(B°). To a good approximation, the SM 
predicts Sf = — £j sin2<£i and Af — for both b — > ccs 
and b — > sqq transitions, where £f = +1( — 1) corre- 
sponds to CP-even (-odd) final states. Recent theoret- 
ical studies within the SM framework [3] find that the 
effective sin20i values, sm2<fi\ s , obtained from b — > 
sqq are expected to agree within O(0. 01) with sin20i 
from the b — ► ccs transition. A comparison of CP- 
violation parameters between these theoretically clean 
b — > sqq modes and b — > ccs decays is an important 
test of the SM. On the other hand, Sf is expected to 
be small for b — > dqq transition because the weak and 
quark mixing phases cancel [3]. In this paper, we report 
measurement of time-dependent CP-asymmetries of 
penguin dominated B meson decays. Among the final 
states, <j>K%, r}'K%, uK° s , K°ir° and J/^K° S are CP 
eigenstates with = —1, while 4>Kj\, rfK®, foKg, 
K%ir Q ir Q , K° S K Q S K% and J/ipK% are CP eigenstates 
with ^ = +1. Since P° -> K + K~K% is a CP-even 
and -odd mixture, of both it has — — (2/+ — 1), 
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where /+ is the CP-even fraction, measured to be 
0.93 ± 0.09 ± 0.05 assuming isospin relation [5]. The 
CP asymmetry parameters for b — ► dqq transition are 
measured from B° — > KgKg decay mode. 

2 Measurement of CP violation 
parameters 

At the KEKB energy-asymmetric e + e _ (3.5 on 8.0 
GeV) collider [5], the T(45) is produced with a Lorentz 
boost of /?7 = 0.425 nearly along the electron beam 
line, (z-axis). Since the B° meson pair is approximately 
at rest in the T(4S) center-of-mass system (cms), At 
can be determined from the displacement in z be- 
tween the fcp and /t ag decay vertices :At ~ (zcp — 
^ta g )/(/37c) ee Az/((3jc). The Belle detector [7] is a 
large-solid-angle magnetic spectrometer that consists 
of a silicon vertex detector (SVD), a 50-layer central 
drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement 
of time-of-fiight scintillation counters (TOF), and an 
electromagnetic calorimeter comprised of CsI(Tl) crys- 
tals (ECL) located inside a super-conducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to 
detect K\ mesons and to identify muons (KLM). 

The intermediate meson states are reconstructed 
from the following decays: 7r° — > 77, Kg — » tt + it~ (de- 
noted by Kg~ hereafter) or 7r°7r° (denoted by K s ° 
hereafter), r\ — > 77 or 7r + 7r~7r°, p — > 7r + 7r~, 77' — » 
p°7 or r]ir + ir~, u — > 7T + 7r~7r°, fo — > 7r + 7r~, <p — > 
K+K- and J/ijj -► £+r(£ = /j,e). We use all com- 
binations of the intermediate states except for the fol- 
lowing cases: r\ — > 7r + 7r~7r° candidates are not used for 
B Q — > rj'Kg decays; 77' — > p°7 candidates are not used 
for P° — > decays; .Kg candidates are not used 

for B° -> J/^iCg, /o-K"s, w^g, ^tt , ^57r 7r and 
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K Q s Kl decays. We reconstruct the B° -> K° s K%Kl 
decays in the K^~K^~K^~ or Kj~ K+~ Kf final 
states. In addition, <j> — » decays are used for 

the B° -> (f>K+~ mode. 

For reconstructed B° — ► /cp candidates without 
a meson, B meson decays are identified using the 
energy difference AE = E™ s — E^ m and the beam- 
energy constrained mass M^ c = \/{E^ m ) 2 — (p™ ls ) 2 , 
where E™ s and p c ^ s are the cms energy and momen- 
tum of the reconstructed B candidate, respectively. 
The signal candidates are selected by requiring 5.27 
GeV/c 2 < Mbc < 5.29 GeV/c 2 and a mode-dependent 
AE window. Only M^ c is used to identify the decay 
B° -> <j)K% followed by <j> -► Other candi- 

date B° — » /cp decays with are selected using 
pf 18 , requiring 0.2 GeV/c < pg 118 < 0.45 GeV/c for 
B° -► J/i/'^l candidates and 0.2 GeV/c < p^ 8 < 0.5 
GeV/c for the others. 

The dominant background for this analysis comes 
from continuum events: e + e~ — > uu,dd,ss, or cc. To 
distinguish these topologically jet-like events from the 
spherical-like £> decay signal events, we combine a set 
of variables [5JH] that characterize the event topology 
into a signal (background) likelihood variable £ s i g (bkg) > 
and impose loose mode-dependent requirements on the 
likelihood ratio TZ s / h = £ s i g / (£ s i g + £bk g )- 

The contributions from BB events to the back- 
ground for B° — > fcp candidates with a K\ are esti- 
mated with Monte Carlo (MC) simulated events. The 
small BB contribution to the background in the B° — > 
rj'Kg mode is also estimated using MC events. We re- 
ject KgKgKg candidates if they are consistent with 
B° -» XcoKs -> {K S K° S )K° S or B° D°K° S -» 
(KgKg)Kg decays, i.e., if one of the KgKg pairs is 
consistent with the Xco mass or D° mass. We estimate 
the contribution in the B° — » 4>Kg sample from B° — > 
K+K'K% and B° -> / #s(/o -> K+R-) decays 
from the Dalitz plot for B° — > K + K~K candidates 
with a method that is described elsewhere [10]. The 
fraction of B° -> K + K~K S events in the i?° -> </>if° 
sample is 2.75±0.14%. The background from i?° — > 
f K s decay for 5° -► K+ K~ K% and 5° -► 4>K° is 
found to be consistent with zero and the influence is 
treated as a source of systematic uncertainty. 

The 6-flavor of the accompanying B meson is iden- 
tified from inclusive properties of particles that are not 
associated with the reconstructed B° — > fcp decay. 
We use two parameters, the 6-flavor charge q and r, to 
represent the tagging information [llj . The parameter 
r is an event-by-event, MC-determined flavor-tagging 
dilution factor that ranges from r = for no flavor 
discrimination to r = 1 for unambiguous flavor assign- 
ment. 

The vertex position for the fcp decay is recon- 
structed using charged tracks that have enough SVD 
hits [12] . The /t ag vertex is obtained with well recon- 
structed track that are not assigned to fcp- A con- 
straint on the interaction-region profile in the plane 
perpendicular to the beam axis is also used with the 
selected tracks. 



We determine CP parameters Sf and Af for each 
mode by performing an unbinned maximum-likelihood 
fit to the observed At distribution. The probability 
density function (PDF) expected for the signal dis- 
tribution, V S ig(At;Sf,Af,q,Wi,Awi), is given by the 
time-dependent decay rate of signal incorporating the 
effect of incorrect flavor assignment. The distribution 
is convolved with the proper-time interval resolution 
function lZ S i g (At), which takes into account the fi- 
nite vertex resolution. The resolution and wrong-tag 
fractions are determined by a multi-parameter fit to 
the At distribution of a high-statistics control sam- 
ple of semileptonic and hadronic b — > c decays [511TB]. 
We determine the following likelihood for each event: 
Pi = (1 - /oi) / [f sig V sie (At')R sig (AU - At 1 ) + (1 - 

/ sig )nkg(^o^bkg(^-^')]rf(^')+/oi^oi(^).Thc 

signal probability / s i g depends on the r region and is 
calculated on an event-by-event basis as a function 
of the following variables: AE and M^ c for B° — > 
J/ipK° s ; p™ 8 for B° -> J/ipK'l and <f)Kj\] M bc and 
K s/h for B° -> 0(-> K° S K° L ) K o. AEj Mbc and Ks/h 
for the other modes. Vhkg(At) is a PDF for background 
events, which is convolved with the background res- 
olution function i?bkg- The term P \(At) is a broad 
Gaussian function that represents a small outlier com- 
ponent [8lll3|. The Sf and Af are determined by max- 
imizing the likelihood function L — JT- Pi(Ati; Sf, Af) 
where the product is over all events. Table [JJ summa- 
rizes the fit results for sin 2(j>1 and A / . Figures [TJ and 
[2] show the At distributions and asymmetries for good 
tag quality events of r > 0.5. The dominant sources of 
systematic error for sin 2(/>f ff stem from the uncertain- 
ties in the resolution function and in the background 
fraction. The dominant sources for Af are the effects of 
tag-side interference (TSI) [14] , the uncertainties in the 
background fraction, in the vertex reconstruction and 
in the resolution function. We study the possible cor- 
relations between 7\L s /b, p™ 8 and r PDFs used for §K\ 
and T]'K^, which are neglected in the nominal result, 
and include their effect in the systematic uncertain- 
ties in the background fraction. Other contributions 
come from uncertainties in wrong tag fractions, the 
background At distribution, rgo and Am e i- A possible 
fit bias is examined by fitting a large number of MC 
events and is found to be small. The dominant sources 
of systematic errors for the B° — > J/tpK° mode are 
the uncertainties in the vertex reconstruction, in the 
resolution function, in the background fraction, in the 
flavor tagging, a possible fit bias, and the effect of the 
TSI. Other contributions are negligible. We add each 
contribution in quadrature to obtain the total system- 
atic uncertainty. The systematic errors are summa- 
rized in Table [2 



3 Summary 

For the B° — > r/'K mode, we determine the statistical 
significance from the obtained statistical uncertainties, 
taking into account the effect of the systematic un- 
certainties. The Feldman-Cousins frequentist approach 
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Table 1. Number of signal -/V s i gna i and results of the fits to At distributions. The first errors are statistical and the 
second errors are systematic. The third error for sin 2<j)f of K+K'Kl mode is an additional systematic error arising 
from the uncertainty of the £f = +1 fraction. 



Mode 


signal 


sin2^f ff 


A S 


<f>K° 


307±21 (Kg), 114±17 (Kl) 


+0.50 ±0.21 ±0.06 


±0.07 ±0.15 ±0.05 


n'K° 


1421±46 (K° s ), 454±39 (K°) 


±0.64 ±0.10 ±0.04 


-0.01 ±0.07 ±0.05 


ujK% 


118±17 


±0.11 ±0.46 ±0.06 


-0.09 ±0.29 ±0.06 


K° s n° 




±0.33 ±0.35 ±0.08 


-0.05 ±0.14 ±0.05 


K S 7r°n° 


307±32 


-0.43 ±0.49 ±0.09 


-0.17 ±0.24 ±0.05 


foK° s 


377±25 


±0.18 ±0.23 ±0.11 


-0.15 ±0.15 ±0.07 


K%K%K% 


185±17 


±0.30 ±0.32 ±0.08 


±0.31 ±0.20 ±0.07 


K + K~K° 


840±34 


±0.68 ± 0.15 ±0.03tg;?ij 


-0.09 ±0.10 ±0.05 


K%K% 


58±11 


S f = -0.38 ± 0.77 ±0.08 


-0.38 ±0.38 ±0.05 


J/4>K° 


7484±87 (K° s ), 6512±123 (K° L ) 


sin 20i = ±0.642 ± 0.031 ± 0.017 


±0.018 ±0.021 ±0.014 




Table 2. Dominant source of systematic error. 



<f>K° V 'K° loK° s K° S 7T° K° S 7v tv° f°K° s K° S K° S K° S K+ K~ K° K° S K° S J/4<K< 

sin 2(f>f 

Resolution function 0.04 0.04 0.05 0.07 0.04 0.02 0.05 0.08 0.06 0.006 

Background fraction 0.04 0.02 0.04 < 0.01 0.05 0.04 0.06 0.01 0.04 0.006 

Background At shape 0.01 < 0.01 < 0.01 0.05 0.05 0.09 0.01 < 0.01 0.04 0.001 



As 



Resolution function 


0.02 


0.01 


0.02 


< 0.01 


0.02 


< 0.01 


0.02 


0.05 


< 0.01 


0.001 


Background fraction 


0.04 


0.02 


0.02 


< 0.01 


0.03 


0.03 


0.06 


0.07 


0.02 


0.009 


Vertex reconstruction 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.009 


TSI QI] 


0.03 


0.02 


0.04 


0.04 


0.04 


0.04 


0.04 


0.03 


0.03 


0.009 
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Fig. 2. Asymmetries for event with good tag (r > 0.5) for B° -> f K%, K + K~K°, K° s -k° , u)K° s , ifgvrV and K%K%. 
Dashed lines show the SM expectation from B° — > J/ipK° measurement. 



[T5] gives the significance of CP violation that is equiv- 
alent to 5.6 standard deviations for a Gaussian er- 
ror. We conclude that we have observed mixing-induced 
CP violation in the mode B° — > rj'K . We do not find 
any significant difference between the results for each 
individual b — > sqq and b — > dqq mode and those pre- 
dicted from SM. Since some models of new physics pre- 
dict such effects, our results can be used to constrain 
these models. However, many models predict smaller 
deviations that we cannot rule out with the current 
experimental sensitivity. Therefore, further measure- 
ments with much larger data samples are required in 
order to search for new, beyond the SM, CP-violating 
phases in the b — » s transition. 
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